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I INTRODUCTION 


BACKGROUND 

This study of the water quality of the Oconee River Basin was conducted 
as part of the expanded scope Flood Plain Information (XFPI) study done 
by the Savannah District Army Corns of Engineers. The Savannah District 
undertook this first pilot XFPI study In which qeographic data banks 
were used as the basis for simulating watershed hvdrology, and the com¬ 
putation of expected annual flood damages. Environmental considerations 
were also Included in the original pilot study objectives but mainly for 
appraisal of wildlife habitat and tradeoffs between the desirability of 
certain land uses. The primary Intent of the XFPI analysis was to analyze 
the effects of alternative futures of the Oconee River Basin development 
In a systematic manner such that realistic comparisons of flood hazards, 
flood damages, and environmental quality could be made between existing 
and alternative future watershed develownent patterns. 

After the viability of the XFPI geographic data bank methodology had been 
successfully demonstrated for flood hazard and damage computations, the 
Savannah District requested the HEC to perform an Oconee River Basin 
water quality study consistent with the XFPI objectives and methodology. 
Thus, the existing and alternative future water quality of the Oconee 
River, within the study area, would be simulated using the geographic 
data bank as the basis for land use inputs to existing HEC water quality 
simulation models. The Storaqe, Treatment, Overflow, Runoff Model (STORM) 
[1] would be used for determining the quantity and quality of land surface 


1 





runoff and dry weather flow and the Water Quality for River/Reservoir 
Systems Model (WORRS) [2] would be used to simulate water quality In the 
river network. That is, the land surface runoff from STORM would be Input 
to WQRRS which would combine all Inflows to the Oconee River and simulate 
the resultant river water quality. 

Some historical data about the water quality of the Oconee River Basin 
were said to be available but the extent and appropriateness of those data 
for these modeling purposes were not known. The HEC undertook the water 
quality studies expecting that the historical data would be satisfactory 
for model calibration for existing conditions. If some aspects of these 
data were not sufficient, then general experience from other water quality 
studies would be used to ascertain acceptable performance of the simulation 
models. 

The Savannah District also requested the HEC to study the land surface 
erosion simulation possibilities using the geographic data bank. In 
particular. It was desirable to erode and transport sediment on a grid 
cell basis as defined by the topography stored In the geographic data bank. 
A new computer program would be developed to Implement this proposed 
methodology. 


SCOPE AND OBJECTIVES 

The objective of the HEC study was to Investigate the applicability of the 
HEC water quality simulation models, STORM and WQRRS, for usage in XFPI 
studies. This was to be accomplished through an evaluation of the water 
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quality Impacts of existing and alternative future land use development 
patterns In the Oconee River Basin. The new HEC grid cell sediment transport 
model would be evaluated In a similar manner. The methodologies for analysis 
of water quality and sediment transport were to be consistent with the 
philosophy (l.e., use of geographic data banks) of the ongoing XFPI pilot 
study. 

The land use data required by the STORM model were to be obtained from the 
geographic data bank through the Hydrologic Parameters (HYDPAR) utility 
program. STORM would also access the land use data for computation of 
sanitary sewage flows. Certain changes were reouired In STORM to utilize 
the grid cell data bank. The HEC would calibrate the models on existing 
data and use the calibrated models to simulate alternative future watershed 
developments. 

STUDY TEAM 

This study was to be carried out entirely by the HEC with minimal direct 
Involvement of Savannah District personnel. The district was to provide 
general guidance about the XFPI study, the objectives of this study, and 
supply existing water quality and sediment data to the HEC. Because this 
study was a special investigation of the applicability of new techniques 
in support of the XFPI study, the district did not feel it was necessary 
to have this water quality modeling expertise developed within their staff. 
The district staff could be trained in the use of this methodology at a 
later date if warranted by their needs and the results of this study. 
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The HEC conducted the study as a team effort with Mr. Jess Abbott being 
responsible for the application of the STORM model. Dr. Michael Gee being 
responsible for the grid cell sediment transport Investigations, and Mr. 

R. G. Willey being responsible for the application of the WORRS model. 
Messrs. Darryl Davis and Pat Webb provided guidance on the XFPI methodology 
and the utilization of the geographic data bank. Mrs. Marilyn Hurst and 
Mr. Paul Ely performed most of the detailed tasks involved to complete the 
project. The drafting was done by Mr. Roger Nutter. 
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II SUMMARY AND CONCLUSIONS 


SUMMARY 


The study objectives were carried out as proposed. The development, 
modification, and implementation of the mathematical models provided the 
means to simulate exlstlnq and assess the future water quality and sediment 
transport characteristics of the Oconee River Basin study area. The STORM 
and WQRRS water quality modeling methodology was successfully Implemented; 
however, the lack of adequate data to calibrate these models made their 
application difficult and the results unsubstantiated except for general 
comparisons with experience from similar studies. 


The grid cell land surface erosion and sediment transport methodology was 
not entirely successful because of problems In using the topography file 
of the geographic data bank. The erosion/sediment transport methodology 
was based on continuously downhill sloping grid cells to the stream 
collection network. Upon application of the method using the Oconee 
topography file, it was learned that the grid cells did not slope con¬ 
tinuously downhill to the collection channels. A concentrated effort was 
made to incorporate the required slope continuity into the geographic data 
bank. This problem could not be resolved within the scope of this 
investigational effort. If further development of this methodology Is 
desired, the erosion/sediment transport methodology must be modified to 
accept the existing data bank topography or the topographic data must 
be edited to conform to the slope-continuity assumptions made In the 
mathematical model. 











The STORM and WQRRS Interfaces with the data bank were developed 
simultaneously with assembly of historical data required to operate and 
calibrate the models. It soon became evident that very few historical 
water quality measurements had been made in the Oconee River and its 
tributaries. Thus, the calibration of the models would have to be based on 
theory and experience. Some water auality measurements had been made and 
these data were used as much as possible. No data were available for the 
important storm runoff periods durlnq which land surface pollutant washoff 
occurs. Concentrations of pesticides, heavy metals and other parameters were 
not specifically mentioned and not evaluated in this study. Such parameters 
cannot be simulated presently by STORM and WQRRS. 

The decision was made to continue with the STORM-WORRS modeling effort with¬ 
out adequate data. This was done mainly so that the general methodology 
could be developed. Had it not been for the desire to develop-and demonstrate 
the STORM-WQRRS methodology for XFPI studies, a more simple water quality 
study would have been recommended for the Oconee study. The recoirmended 
water quality study methodology would have been commensurate with the 
detail of existing data and the degree of detail warranted by the study 
objectives. 

Assumptions were made about the basic water quality Inputs from the Upper 
Oconee River (outside the study area) and about the land surface runoff 
within the study area. Because there were no data on the river water quality 
during storm events, the calibration of reasonable values of land surface 
runoff and Incoming river water quality required much more time than 
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anticipated. The land surface runoff simulation was reviewed In some detail 
for Reaches 1 and 2, Figure III-3. Upon achieving reasonable results for these 
reaches. It was assumed that subsequent reaches simulated in a similar manner 
would also have reasonable results. Thus, the other basins were simulated 
and it was not until all of the basins were aggregated in the WQRRS receiving 
water model that it was apoarent some basin results were unreasonable. The 
complexity of the WORRS model did not facilitate the timely appraisal of these 
potential problems. At that late date in the study, some of the basin land 
surface runoff simulations had to be rerun and the river system was simulated 
again. The new results were acceptable. The land surface runoff from the 
basins should have been reviewed more thoroughly during the initial simulations; 
however, such review was not emphasized because of the data limitations. 

The existing land use condition was simulated using non-point source land 
surface runoff, point sources within the tributary subbasins, and point source 
Inputs from the two main sewage treatment plants in the study area. The 
treatment plant loads appear to cause the most significant Impact on the 
North Oconee River for existing conditions. The pollutant loadings from 
the sewage treatment plants on the North Oconee contribute approximately 
80-93% (depending on the parameter) of the total loads (point and non-point 
sources). Therefore further reduction in the loads from sewaqe effluent 
would appear to have the greatest effect in upgrading the water quality for 
existing conditions. There does not. appear to be a significant water 
quality problem under existing conditions because, with the possible ex¬ 
ception of conforms, established water quality standards are not exceeded. 










The basins having the most significant Impact are 1A, 2, and 3 on the North 
Oconee and 6A, 6B, and 16 on the Middle Oconee, Figures III-l and II1-2. 

The sewage treatment plant effluents on the Middle Oconee contribute 48-80% 
of the total loads reaching the Middle Oconee above the confluence. The 
concentrations do not exceed established water quality standards under the 
existing development with the exception of coliforms. The most irinedlate 
Improvements could be made by reducing the loads from treatment plants. 

The North Oconee River watershed contributes a significantly larger pollutant 
load to the Main Oconee River than does the Middle Oconee River. The North 
Oconee watershed has significantly larger loadings from the sewage treatment 
plants. Specifically the North Oconee plant contributes about 70% of the 
total load (point and non-point sources) to the Main Oconee River. In 
addition, effects of the loadings from the treatment Diants are much more 
pronounced In the North Oconee than the Middle because the natural flows 
from the upstream watershed on the North Oconee provide much less dilution. 

The future water quality for the Oconee River study area was simulated for 
the Alternative C, Table III-l, land use development plan. In general, the 
sources of water degradation are the same as those defined for existing 
conditions. On the North Oconee the contributions due to treatment plant 
loadings ranged from 81-88% of the total loadings and on the Middle Oconee 
the plant loadings contribute 48-88% of the totals. 

i 

j 

While the percentage contribution from sewage treatment loads appears to ] 

remain constant, the total loads have increased somewhat from existing to 
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the alternative future C. The increases are not major due to the relatively 
small percentage change in land use change in the total study area. The 
major Impacts were shown to be In subbasins 1A, 2, and 3 on the North Oconee 
and 6A, 6B, and 16 on the Middle Oconee because these subbasins experienced 
the greatest degree of urbanization. These subbasins and the loadings from 
the 2 main treatment plants tend to create "shock loadings" In the reaches 
immediately downstream of the effluent outfalls or tributary Inflows. 
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CONCLUSIONS 


Conclusions with respect to both the technical feasibility and the 
suitability of the methodology for water quality and land surface erosion/ 
sediment transport studies in support of the XFPI program will be made in 
this section. The impact of the future land use plans on Oconee River 
water quality were discussed in the preceding Summary Section. The study 
objectives, budget, and availability of data to support the proposed study 
methodology are important factors in determining the most appropriate 
technology for a project. The STORM and WQRRS water quality simulation 
models have been shown to be technically feasible for water oualitv studies. 
The appropriateness of their use for XFPI studies is a much more important 
question which can only be determined by the water quality objectives of 
each study application. The STORM and WORRS computer programs and attendant 
study methodoloqy are quite appropriate for detailed water quality studies. 
That is, these models provide a good simulation of the ohysical water 
quality svstem, both for land surface runoff and receiving waters. The 
STORM model provides a relatively simple simulation of land surface runoff. 
The WQRRS model performs a rather complex simulation of receiving water 
quality and requires much more comprehensive input data than does STORM. 

The general objective of XFPI studies Is the analysis of the hvdroloqic, 
economic, and environmental impact of future land use development patterns. 
To accomplish this, the existing system must be reoresented satisfactorily 
in simulation models. There must also be a consistent, logical means to 
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generate and compare alternative futures. Because there are so many unknowns 
with respect to specific location and type of future land use patterns, river 
regulation, and waste water management facilities, the analysis of futures 
can be less detailed than known conditions. Methods should capture the 
esseoce of the future conditions without being overly complex about the 
specific types and locations of the development. 

In accordance with the above objectives, the STORM model provides both the 
type of Information and the level of technical detail which are appropriate 
for XFPI studies. The basic land use parameters of the STORM model are 
readily derived from the geographic data bank. Other Input to the STORM 
model can be easily obtained or estimated from previous experience. 

For application of the STORM model, continuous rainfall and runoff data are 
recommended for a period of several years to calibrate the hydrologic parameters 
of the model. Pollutant loadings at the subbasin outlets should be measured 
throughout several major storm events during the multiyear hydrologic calibration 
period. If these data are not available, the STORM model should not be used 
unless acceptable results can be obtained from use of coefficients derived from 
similar studies. The availability of data should be determined early in the 
study so that data collection efforts cannbe arranged as necessary. 

The water quality simulation capability of the WQRRS model seems to be more 
comprehensive than required for the general XFPI study. The many data 
requirements limit the utility of this model. The WQRRS model would be 
required if a more comprehensive understanding of the water quality condition 
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Is desired. This might be the case for reservoir regulation studies or 
major river/reservoir studies with specific water quality objectives. For 
the XFPI level of complexity, a more simple receiving water analysis seems 
appropriate. An analysis conmensurate with the complexity of the STORM 
model Is desirable. The HEC Is presently developing that type of simplified 
receiving water model. 

If a detailed water quality study had been required for the Oconee XFPI study, 
then data collection efforts should have been started as soon as It was 
determined that the historical data were Inadequate. The collection and 
analysis of field data would have required a considerably larger study budqet, 
on the order of 8 to 10 times the Initial water Quality study budget of 
$ 20 , 000 . 

For application of the WORRS model, climatic, pollutant Dolnt-source loading 
data, and/or results from STORM must be known for the entire calibration 
period. Most Importantly, in-stream water quality measurements must be 
available for several major parameters both during storm runoff periods 
(preferably the same periods as used for STORM) and low flow or other critical 
water quality periods. These data should be available for at least one 
location in the river system depending upon how much variation there Is In 
the land use and stream regimes in the river network. If only one location 
were available. It should be at the downstream boundary of the river network 
so that the Integrated effects of the land surface runoff and In-stream 
quality changes are measured. 
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Ill OCONEE RIVER SYSTEM DATA AVAILABILI TY 
GENERAL 

The Oconee River begins in the Georgia counties of Barrow and Jackson north 
of Athens and flows south through the middle of Georgia. After It joins the 
Ocmulgee River near Hazlehurst, it becomes the Altamaha River which flows 
southeast to the Atlantic Ocean. A location map is shown in Figure III-1. 

The study boundaries for this project include the Oconee River drainage be¬ 
tween the Currey Creek dam site on the North Oconee and the State Highway 
Bridge 33 on the Middle Oconee down to a location 8 miles below Barnett 
Shoals Dam on the Main Oconee (i.e., inflow to Wallace Reservoir). The 
Oconee River study area is shown in Figure 111-2 and schematics of the study 
area are shown in Figures 111-3 and III-4. The schematics include the 
location of all major tributaries, sewage treatment plant effluents, and 
the Athens water supply intake. 

The historical period to be used for analysis was selected using the 
following criteria: 

(1) A low flow period. 

(2) A period with several significant rain events. 

(3) A relatively recent period (i.e., existing conditions). 

(4) A short duration (i.e., one month). 

Water quality sampling points were located on the Middle Oconee River in 
1970 and on the North Oconee in 1974. During 1970-1975, the dryest one 
month period having 3 to 4 significant rain events was October 1970. This 
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period has as much available water quality data as any period in 1970 to 
1975 and was therefore selected for analysis. 

The various aspects of specific data availability will be discussed in the 
remainder of this chapter. 

























NORTH OCONEE 
BOUNDARY INFLOW 
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Sub basin 12C (19.0) 




Fig HI-3. Schernatic of North Oconee River 
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MIDDLE OCONEE 
BOUNDARY INFLOW 
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Subbasin 5C (31.2) 



Subbosin 5D —Bear Cr. (29.5) 

Subbosln 5A— Turkey Cr.(25.9) 

Subbosin 5B(24.6) 

Middle Oconee USGS Gage 2-2175(22.3] 
Subbosin 6B(22.Q) 

Subbosin 6A (19.9) 

Middle Oconee STP(17.3) 
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Middle and Main Oconee River 
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Subbasin 15 (9.8) 



Fig. m-4. Schematic of Middle and Main Oconee River 




























METEOROLOGY 


The weather data for the analysis were obtained from the National Weather 
Service using the Athens Municipal Airport Weather Station. A magnetic 
tape of hourly rainfall at five stations for the STORM model input and 
another tape of dry and wet bulb air temperature, barometric pressure, 
wind speed and cloud cover for the WORRS model input were obtained from 
the Ashville, North Carolina office of the U.S. Weather Service. 

These input data were the easiest to obtain of all the required input 
data for either STORM or WQRRS. See Chapter IV for discussion of models. 
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LAND USE 


Land use Is one of the most Important input variables for STORM. It is 
especially important in this study since one of the main objectives is to 
assess the impact of future development (as characterized by land use) on 
the water quality of the Oconee River. Land use for each STORM watershed 
was taken directly from output from HYDPAR (an HEC-developed utility program 
to calculate hydrologic parameters from a grid cell data ban!;). The specific 
land use categories that were used in this study are as follows: 


Code No. 
1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


Designation 

Developed Open Space (lawns, parks, golf courses, 
cemeteries and rights-of-way) 

Low Density Residential 

Medium Density Residential 

High Density Residential 

Agriculture (cultivated land, row crops, small grain) 
Industrial 

Commercial (strip and isolated commercial) 

Pasture 

High Density Commercial (downtown areas and shopping centers) 

Institutional 

Natural 


Some alterations were made to a few of the land use categories. Low 
Density Residential and Medium Density Residential were combined into a 
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single category. Hardwoods, pines, and wetlands were combined Into a single 
category (Natural). Roads, land fills and water bodies were not simulated, 
however, the areas of each of these categories were subtracted from each 
watershed area. Table III-l shows the land use for each STORM watershed 
for both existing and one alternative future (1990C). 














TABLE m-1 


LAND USE BY WATERSHED AND REACH 


REACH 

NO. 

SUB-BAS IN 
NO. 

AREA 

(ACRES) 

LAND USE 

EXIST 

* 

k> 

mm 

1 

20 

1870 







Low Residential 

1.1 

1.1 




High Residential 

0.0 

0.0 




Commercial 

0.0 

0.0 




Industrial 

0.0 

0.0 




Pasture 

25.5 

25.5 




Natural 

73.4 

73.4 

1 

17 

3700 







Low Residential 

2.3 

2.3 




High Residential 

0.0 

0.0 




Conmerclal 

0.0 

0.0 




Industrial 

0.0 

0.0 




Pasture 

38.4 

38.4 




Natural 

59.3 

59.3 

I 

18 

5620 







Low Residential 

1.8 

1.8 




High Residential 

0.0 

0.0 




Commercial 

0.0 

0.0 




Industrial 

0.0 

0.0 




Agricultural 

0.8 

0.8 




Pasture 

52.4 

52.4 




Natural 

45.0 

45.0 

1 

19 

2180 







Low Residential 

1.1 

1.1 




Medium Residential 

2.4 

2.4 




Commercial 

0.0 

0.0 




Industrial 

0.0 

0.0 




Pasture 

33.9 

33.9 




Natural 

62.6 

62.6 

1 

IB 

3741 







Low-Medium Residential 

2.3 

6.2 




High Residential 

0.4 

0.4 




Conmerclal 

0.3 

0.6 




Industrial 

0.0 

0.0 




Agricultural 

31.3 

28.2 




Pasture 

7.5 

7.4 




Natural 

57.8 

54.8 




Open 

0.4 

0.6 


NOTE: Water Area Is Included In Natural 
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TABLE III-l 


LAND USE BY WATERSHED AND REACH 


REACH 

NO. 

SUB-BASIN 

NO. 

AREA 

(ACRES) 

LAND USE 

EXIST 

l 


1 

1A 

2928 







Low-Medium Residential 

11.3 

13.0 




Hlqh Residential 

0.6 

1.8 




Commercial 

0.2 

0.6 




Industrial 

4.9 

11.2 




Agricultural 

9.7 

6.3 




Pasture 

11.1 

7.9 




Institutional 

0.2 

2.8 




Natural 

60.0 

51.3 




Open 

2.0 

5.0 

1 

2 

41254 







Low-Medium Residential 

2.1 

3.1 




High Residential 

0.1 

0.2 




Commercial 

0.1 

0.2 




Industrial 

0.0 

0.0 




Pasture 

41.8 

41.7 




Institutional 

0.1 

0.1 




Natural 

50.9 

50.0 




Open 

0.3 

0.4 

2 

3 

2272 







Low-Medium Residential 

30.9 

33.1 




High Residential 

5.9 

4.4 




Commercial 

4.3 

4.6 




Industrial 

14.4 

15.1 




Agricultural 

2.7 

2.1 




Pasture 

3.1 

2.3 




Hlqh Commercial 

1.0 

1 .0 




Institutional 

1.6 

1.5 




Natural 

28.8 

26.6 




Open 

7.3 

6.7 




Roads 

0.0 

2.6 
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TABLE III-l 


LAND USE BY WATERSHED AND REACH 


REACH 

NO. 

SUB-BASIN 

NO. 

AREA 

(ACRES) 

LAND USE 

EXIST 

W 

to 

1990C 

% 

2 

4 

7915 







Low-Medium Residential 

20.0 

21.5 




High Residential 

2.6 

2.6 




Corrmercial 

1.8 

2.0 




Industrial 

1.5 

15.3 




Agricultural 

24.0 

17.9 




Pasture 

6.3 

4.1 




Hiqh Commercial 

0.0 

0.5 




Institutional 

0.7 

1.0 




Natural 

43.1 

31.8 




Open 

0.0 

0.6 




Roads 

0.0 

2.9 

2 

12B 

1756 







Low-Medium Residential 

18.6 

18.6 




Hiqh Residential 

4.q 

6.0 




Commercial 

4.2 

3.8 




Industrial 

o.n 

0.0 




Agricultural 

12.3 

9.7 




Pasture 

0.6 

0.3 




High Commercial 

3.5 

3.5 




Institutional 

34.6 

34.1 




Natural 

1Q.9 

18.7 




Open 

1.4 

2.7 




Roads 

9.0 

2.7 
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TABLE III-l 


LAND USE BY WATERSHED AND REACH 


REACH 

NO. 

SUB-BASIN 

NO. 

AREA 

(ACRES) 

LAND USE 

EXIST 

? 

1990C 

l 

2 

12C 

2702 







Low-Medium Residential 

17.7 

19.7 




High Residential 

3.2 

4.6 




Commercial 

3.8 

3.7 




Industrial 

2.4 

15.2 




Agricultural 

8.2 

5.3 




Pasture 

10.0 

5.5 




High Commercial 

1.1 

0.9 




Institutional 

3.2 

2.8 




Natural 

39.0 

31.5 




Open 

11.4 

7.2 




Roads 

0.0 

3.5 

2 

12A 

3191 







Low-Medium Residential 

7.9 

9.0 




High Residential 

0.5 

2.8 




Commercial 

1.3 

1.3 




Industrial 

0.3 

1.9 




Agricultural 

11.8 

9.0 




Pasture 

5.1 

3.7 




High Commercial 

0.0 

0.7 




Institutional 

3.4 

5.8 




Natural 

65.9 

60.5 




Open 

3.8 

4.8 




Roads 

0.0 

0.4 

3 

5C 

3456 







Low Residential 

1.6 

1.5 




High Residential 

0.4 

0.4 




Commercial 

0.5 

0.7 




Industrial 

0.0 

0.0 




Agricultural 

16.6 

16.5 




Pasture 

4.0 

4 .0 




Natural 

75.9 

75.9 




Open 

1.0 

- - - _ 

1.0 


















TABLE III-l 


LAND USE BY WATERSHED AND REACH 



REACH 

SUB-BASIN 

AREA 

LAND USE 

EXIST 

NO. 

NO. 

(ACRES) 


* 

3 

50 

13020 

Low Residential 

2.3 




Medium Residential 

1.2 




Commercial 

0.3 




Industrial 

0.0 




Agricultural 

0.9 




Pasture 

50.3 




Institutional 

0.2 




Natural 

43.9 




Open 

0.9 

3 

5A 

5723 

Low-Medium Residential 

18.3 




High Residential 

1.7 




Commercial 

0.6 




Industrial 

1.1 




Agricultural 

9.5 




Pasture 

3.9 




High Commercial 

0.2 




Institutional 

0.6 




Natural 

62.8 




Open 

1.3 

3 

5B 

3385 

Low-Medium Residential 

8.1 




Hlqh Residential 

0.9 




Commercial 

1.1 




Industrial 

0.0 




Agricultural 

17.5 




Pasture 

4.7 




High Commercial 

0.2 




Natural 

67.4 




Open 

0.1 

3 

6B 

2658 

Low-Medium Residential 

23.9 




High Residential 

5.4 




Commercial 

7.2 




Industrial 

0.0 




Agricultural 

0.9 




Pasture 

3.0 




High Commercial 

1.7 




Institutional 

2.7 




Natural 

53.4 




Open 

1.8 





















TABLE III-l 


LAND U c -E BY WATERSHED AND REACH 





It- 

I? 


• 



f; 
* . 


* 

5 

-4 

» v 

fc 

<: 

♦ 


REACH 

NO. 

SUB-BASIN 

NO. 

AREA 

(ACRES) 

LAND USE 

EXIST 

* 

* 

1990C 

f 

3 

6A 

3443 







Low-Medium Residential 

47.6 

49.5 




High Residential 

2.4 

3.2 




Commercial 

5.9 

5.9 




Industrial 

0.0 

0.0 




Agricultural 

1.9 

0.9 




Pasture 

4.9 

3.7 




High Commercial 

1.7 

1.7 




Institutional 

3.0 

2.8 




Natural 

32.0 

29.5 




Open 

0.6 

1.2 




Roads 

0.0 

1.6 

3 

7 

27410 







Low Residential 

3.0 

2.8 




Medium Residential 

0.8 

1.5 




Commercial 

0.1 

0.2 




Industrial 

0.0 

0.2 




Agricultural 

16.0 

15.4 




Pasture 

37.1 

37.3 




Institutional 

0.1 

0.2 




Natural 

42.8 

42.3 




Open 

0.1 

0.0 

4 

8 

10260 







Low-Medium Residential 

13.1 

15.2 




High Residential 

0.8 

1 .4 




Commercial 

1.1 

1.2 




Industrial 

0.0 

0.0 




Agricultural 

16.6 

14.7 




Pasture 

24.8 

25.5 




High Conmercial 

0.0 

0.1 




Institutional 

0.2 

0.5 




Natural 

43.0 

40.5 




Open 

0.4 

0.7 




Roads 

0.0 

0.2 

















TABLE 111-1 


LAND USE BY WATERSHED AND REACH 


REACH 

NO. 

SUB-BASIN 

NO. 

AREA 

(ACRES) 

LAND USE 

EXIST 

to 

H 

4 

9 

1290 







Low Residential 

0.8 

0.8 




Medium Residential 

8.3 

7.8 




Commercial 

1.2 

1.5 




Industrial 

0.0 

0.0 




High Residential 

0.0 

3.4 




Agricultural 

20.6 

17.8 




Pasture 

5.2 

4 .2 




Institutional 

0.0 

1.8 




Roads 

0.0 

2.5 




Natural 

63.9 

60.3 

4 

10 

5946 







Low Residential 

0.5 

0.5 




Medium Residential 

7.2 

9.7 




Commercial 

0.7 

0.8 




Industrial 

0.4 

0.4 




Agricultural 

36.9 

34.7 




Pasture 

1.7 

1 .7 




High Connercial 

0.0 

0.1 




Institutional 

1.3 

1.6 




Roads 

0.0 

0.4 




Natural 

51.3 

50.1 

4 

11 

3460 







Low Residential 

1.5 

1.5 




Medium Residential 

0.1 

0.6 




Commercial 

0.4 

0.4 




Industrial 

0.8 

2.3 




Agricultural 

10.7 

9.8 




Pasture 

13.2 

12.9 




Institutional 

0.1 

1 .0 




Natural 

73.2 

71.4 

4 

13 

3428 







Low-Medium Residential 

20.6 

32.1 




High Residential 

1.3 

2.7 




Commercial 

0.2 

0.5 




Industrial 

0.0 

2.8 




Agricultural 

12.4 

8.5 




Pasture 

8.3 

7.0 




High Conmerclal 

0.0 

0.2 




Institutional 

1.5 

2.3 




Natural 

55.7 

43.1 




Open 

0.0 

1.0 
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TABLE III-l 



REACH 

SUB-BASIN 

NO. 

NO. 




LAND USE 


Low Residential 

Medium Residential 

Conmerclal 

Industrial 

Agricultural 

Pasture 

Natural 


Low Residential 

Medium Residential 

Commercial 

Industrial 

Agricultural 

Pasture 

Institutional 

Natural 


EXIST 



11259 


Low-Medium Residential 

High Residential 

Commercial 

Industrial 

Agricultural 

Pasture 

Institutional 

Natural 

Open 


Low Residential 

High Residential 

Commercial 

Industrial 

Pasture 

Natural 


Low Residential 

High Residential 

Commercial 

Industrial 

Agricultural 

Pasture 

Natural 


0.0 

3.6 

30.4 

27.9 

5.7 

5.1 

0.1 

0.9 

56.6 

51.8 

0.4 

0.4 

0.4 

0.4 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

48.7 


50.9 

50.9 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.5 

0.5 

83.1 

83.1 

16.4 

16.4 
















TABLE II1-1 


LAND USE BY WATERSHED AND REACH 


REACH 

SUB-BASIN 

AREA 

LAND USE 

EXIST 

1990C 

NO. 

NO. 

(ACRES) 



% 

5 

22 

39550 

Low Residential 

1.0 

1.0 



l 

| Medium Residential 

0.2 

0.2 




Commercial 

0.0 

0.0 




Industrial 

0.0 

0.0 



1 

Agricultural 

0.1 

0.1 




Pasture 

23.9 

23.9 




Natural 

74.8 

74.8 

5 

25 

9880 

Low Residential 

0.0 

0.0 




High Residential 

0.0 

0.0 


' 


Commercial 

0.0 

0.0 




Industrial 

0.0 

0.0 




Pasture 1 

15.1 

15.1 




Natural 

| 

84.9 

84.9 

5 

24 

1320 

Low Residential 

0.0 

0.0 




High Residential 

0.0 

0.0 



| 

Commercial 

0.0 

0.0 



1 

Industrial 

0.0 

0.0 




Pasture 

86.7 

86.7 




Natural 

13.3 

13.3 















RIVER GEOMETRY 


Cross section data at irreqular Intervals along the entire stream system 
were provided by the Savannah District. The data were provided in a format 
for input to computer program HEC-2, Water Surface Profiles [3]. HEC-2 
output provided Information on energy grade line elevations which is a 
required Input to WQRRS. More importantly, these same cross sections are 
Input to computer program GEDA, Geometric Elements from Cross Section 
Coordinates [4], GEDA provides output of vertically layered geometric 
data (l.e., cross section area, top width, hydraulic radius, composite 
Manning's n, etc.) at regularly spaced nodal points (e.g., one half mile 
apart), as required by WQRRS. 

The preparation of geometric data for the WQRRS model is relatively 
automatic once the basic data of station-elevation coordinate points 
have been obtained either from field surveys or from contour maps. 









HYDROLOGY 


The Oconee basin was found to have little available hydrologic data. 

For the desired watershed modeling purposes, hydrologic data can be 
considered virtually non-existent. This serious lack of data required 
numerous assumptions. The accuracy of these assumptions cannot be 
evaluated except that the results did not seem unreasonable in terms 
of general hydrologic engineering judgment. 

Only one USGS stage gage with hourly flow records was in operation during 
the selected study period, 1570. This gage is located at river mile 22.3 
on the Middle Oconee River. To obtain the inflow across the study 
boundaries on the Middle and North Oconee Rivers (i.e., river mile 31.5 
and 34.8 respectively), the hourly flow rate at the USGS gage was multi¬ 
plied by the ratio of drainage area above the gage to that above each 
boundary. 

Modified Puls routing criteria for the North, Middle and Main Oconee 
Rivers were provided by the Savannah District for selected control 
points. These data were linearly interpolated to obtain criteria at 
each load point (i.e., tributary inflows, sewage treatment plant 
effluents, and withdrawal locations). These data are shown in Table 
111 - 2 . 
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TABLE II1-2 

MODIFIED PULS ROUTING CRITERIA 


NORTH OCONEE RIVER 


Bate-' 

STORAGE 

(AF) 

OUTFLOW 

(cfs) 

RIVER 

MILE 

STORAGE 

(AF) 

OUTFLOW 

(cfs) 

33.5-34.8 

16 

180 

23.2-23.7 

136 

5282 


25 

360 


253 

10600 


112 

1800 

| 


289 

13921 

31.8-33.5 

27 

180 

21.7-23.2 

206 

5282 


42 

360 


340 

10600 


247 

1800 


425 

13921 

31.0-31.8 

14 

l 

180 

20.3-21.7 

65 

5282 


22 

360 


125 

10600 


119 

1800 


163 

13921 

28.7-31.0 

35 

180 

19.9-20.3 

35 

5312 


56 

360 


67 

10900 


352 

1800 


85 

14338 

26.0-28.7 

55 

180 

19.0-19.9 

63 

5312 


81 

360 


123 

10900 


246 

1800 


155 

14338 

24.3-26.0 

27 

180 

15.6-19.0 

440 

5312 


42 

360 


786 

10900 


174 

1800 


995 

14338 

23.8-24.3 

136 

180 

12.3-15.6 

435 

5312 


253 

360 


804 

10900 


289 

1800 


1039 

14338 

23.7-23.8 

30 

5282 





57 

10600 





65 

13921 
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TABLE 111-2 

MODIFIED PULS ROUTING CRITERIA 


MIDDLE AND MAIN OCONEE RIVER 


BH 

STORAGE 

(AF) 

OUTFLOW 

(cfs) 


STORAGE 

(AF) 

OUTFLOW 

(cfs) 

31.2-31.5 

5 

270 

16.5-17.0 

52 

7100 


7 

525 


100 

13500 


17 

2700 


122 

17000 

29.5-31.2 

521 

7100 

15.1-16.5 

407 

7100 


829 

13500 


785 

13500 


964 

17000 


961 

17000 

25.9-29.5 

1148 

7100 

12.9-15.1 

501 

7100 


2397 

13500 


961 

13500 


2937 

17000 


1391 

17000 

24.6-25.9 

145 

7100 

12.3-12.0 

50 

7100 


279 

13500 


95 

13500 


361 

17000 


120 

17000 

22.0-24.6 

380 

7100 

11.4-12.3 

353 

8675 


636 

13500 


654 

19300 


784 

17000 


819 

2 5222 

19.9-22.0 

395 

7100 

9.8-11.4 

564 

8675 


731 

13500 

i 

! 

1173 

19300 


885 

17000 

i 

1462 

25222 

17.3-19.9 

571 

7100 

8.9-9.8 

214 

8675 


1034 

13500 


364 

19300 


1247 

17000 


433 

25222 

17.0-17.3 

52 

7100 

8.0-8.9 

133 

8675 


100 

13500 


206 

19300 


122 

17000 


235 

25222 









TABLE 111-2 

MODIFIED PULS ROUTING CRITERIA 


MIDDLE AND MAIN OCONEE RIVER 



STORAGE 

(AF) 

OUTFLOW 

(cfs) 

RIVER 

MILE 

STORAGE 

(AF) 

OUTFLOW 

(cfs) 

.0 

398 

8675 

0.8-2.8 

359 

8675 


777 

19300 


1031 

19300 


957 

25222 


1249 

25222 

.2 

734 

8675 

0 . 0 - 0.8 

360 

8675 


1318 

19300 


983 

19300 


1630 

25222 


1248 

25222 

.8 

571 

8675 





1440 

19300 





1742 

25222 













WATER QUALITY 


The Oconee basin was found to have little available water quality data of 
practical use to this study. It was originally thought that adequate water 
quality data were available to calibrate the water quality simulation models, 
STORM and WQRRS. This serious lack of data required numerous assumptions which 
could not be verified with field data. The only evaluation that could be made 
was that the results appeared reasonable In the light of other experience. 

Since 1968, water quality data have been collected at the USGS stage gage 
on the Middle Oconee River (l.e., river mile 22.3) and since 1974 at the 
Athens water Intake on the North Oconee River. The data from these two 
sampling locations together with the data from Smith [5], Appendix A, were 
used to estimate base flow quality data for the boundary condition and the 
tributaries. The boundary quality condition was held constant during 
storm events and the tributary Inflow quantity and quality was obtained 
from the STORM model output. 








IV MODELING CONCEPTS APPLIED 


STORM 

The Storage, Treatment, Overflow Runoff Model (STORM) Is a continuous 
simulation model designed to be used in metropolitan master planning 
studies for evaluating storage and treatment capacities required to 
reduce overflows. Pollutograph (pollutant mass-emmisslon rates) loadings 
can also be computed for use in a receiving water assessment model. 

Since STORM is intended for use in planning studies or for screening 
alternatives, some of its analytical techniques are necessarily 
simplified. For example, the two procedures used to compute the quantity 
of runoff are the coefficient method and the United States Soil Conservation 
Service (SCS) method. In the coefficient method, a single land-use weighted 
runoff coefficient is applied to each hour of rainfall excess above depression 
storage to compute runoff. The runoff coefficient is a function of only the 
respective runoff coefficients for the pervious and impervious areas of the 
watershed. Antecedent conditions and rainfall intensity are not taken Into 
account using this method. 

The SCS runoff curve number technique is considered to be conceptually more 
correct than the coefficient method. The SCS curve consists of a nonlinear 
relationship between accumulated rainfall and accumulated runoff. Since 
STORM requires a continuous analysis, a nrocedure has been added that 
computes the curve number for each event based on the number of dry hours 
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since the previous runoff event and prior evapotranspiratlon and percolation. 
Unit hydrographs can be used to transform the surface runoff excesses Into 
basin outflow hydrographs. 


Loads and concentrations for six basic water quality parameters are computed. 
These are suspended and settleable solids, biochemical oxygen demand, total 
nitrogen, total orthophosphate, and total coliforms. Urban and nonurban 
areas may be described by up to 20 land uses. Other features of STORM are 
the capabilities to compute snowfall/snowmelt, dry-weather flow quantity 
and quality, and land surface erosion. 

STORM has a unique advantage of being able to accept discharge hydrographs 
as input for computing the associated wash-off of constituents. 

WQRRS 

The Water Quality for River-Reservoir Systems (WQRRS) model has 
capability for ecologic evaluation of rivers or reservoirs. It is a 
dynamic continuous simulation model. The model consists of three separate 
but integrable modules. These are the reservoir module, the stream 
hydraulic module, and the stream quality module. Since each module is a 
stand-alone proqram, the reservoir, the streamflow routing, or the stream 
water quality module may be executed, analyzed and interpreted independently. 
The three computer proorams may also be integrated into a complete river 
basin water quality analysis. 

The reservoir section of the proqram estimates the water quality condition 
in deep impoundments that can he represented as one-dimensional systems 

39 


i 


J 










In which the Isotherms, or contours of any parameter, are horizontal. 

This approximation Is generally satisfactory in lakes with lonq residence 
times. However, the approximation is less satisfactory in shallow 
Impoundments or those that have a rapid flow-through time. Systems 
that have a rapid flow-through time are often fully mixed and can be 
treated as slowly moving streams using the stream section of the model. 

The stream hydraulic section of the model includes six hydraulic calcu¬ 
lation options. This module Is capable of handling hydraulic behavior 
for both the "gradually varied" steady and unsteady flow regimes. Peak 
flows from storm water runoff or irregular hydropower releases can be 
accurately represented. 

In the stream quality module, the rate of transport of quality parameters 
can be accurately represented and peak pollutant loads into the steady or 
unsteady hydraulic environment can be simulated. The stream portions of 
WQRRS have two automatic interface options for use with the STORM model. 









V MODELING RESULTS 


STORM RUNOFF QUANTITY AND QUALITY 

The approach used to calculate storm runoff was to sub-divide the total 
study area Into a number of Individual watersheds and apply STORM to each. 
Criteria affecting the number of watersheds Include the degree of refinement 
In discrete points along the receiving water body where Individual 
calculations are to be made and manageability of data for the entire 
study. A total of 32 Individual watersheds were Identified and are shown 
on Figure III-2. Several of the watersheds In the Savannah District Data 
Bank were further subdivided so as to provide better definition of quality 
In the urbanized river reaches. 

The first major effort in the STORM application was to assemble and edit 
hourly precipitation data. Five recording rain gages exist In or near 
the study area; the locations are shown In Figure lll-Z. Only the Athens 
gage (Station No. 0435) was a Class I U.S. Weather Service gage. The 
other four gages are supplemental locations and, as a result, the data 
have not been corrected for gage failures. Numerous places existed on 
the tapes where the gage had failed and the accumulated precipitation was 
shown In the first hour of resumption of recording. A special editor 
program was written to locate these gage failures and redistribute the 
precipitation to the hours in which they occurred. The precipitation 
was distributed evenly over the hours of gage failure. 






The continuous precipitation histories (1948-72) were used to assess the average 
annual land surface erosion for both existing and future conditions. The single 
year 1970 was used for storm water simulations since it was the month of 
October 1970 that was studied in the receiving water analysis using WQRRS. 

No data existed with which to calibrate the rainfall-runoff calculations in 


STORM for the Oconee study. The various soil moisture characteristics required 

for the SCS runoff method were estimated. The October percentage runoff for 

several nearby gages served as a guide. The tributary flows, when combined 

and routed to a gage near the downstream study boundary showed fair agreement 

with the observed, however no check could be made on Individual tributary 

flows. The estimated runoff characteristics are shown below: 

Soil Moisture Max Initial 

at Saturation Abstraction Capacity 

LAND USE (SMAX), Inches (DEPR), Inches 

Open 

14.40 

2.88 

Low D. Residential 

11.70 

1.17 

Medium D. Residential 

7.54 

0.11 

Low-Medium D. Residential 

9.62 

0.18 

High D. Residential 

4.28 

0.06 

Agricultural 

7.54 

1.51 

Industrial 

3.33 

0,05 

Commercial 

2.34 

0.04 

Pasture 

10.00 

2.00 

High D. Commercial 

1.24 

0.02 

Institutional 

3.16 

0.05 

Roads 

4.70 

0.07 


Natural 
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No data existed with which to calibrate the runoff quality calculations In 
STORM for the Oconee Study. The various pollutant accumulation rates required 
to regulate the quality were estimated based on data from the literature 
[6,7]. Adjustments were made during calibration so that tributary storm 
water concentrations for existing conditions did not greatly exceed certain 
measured concentrations in the river. Storm water quality calculations were 
not calibrated directly to the river concentration for three reasons: 

1) minimal data existed, 2) the data consisted of grab samples taken at 
Infrequent Intervals, and 3) there were no indications that the measurements 
were taken during periods of tributary storm runoff. Table V-l shows the 
adopted pollutant accumulation rates. 


TABLE V-l 

Pollutant Accumulation Rates (lb/ac/day) 


Land 

Use 

Susp 

Solids 

11 

mum 

B0D 5 

N 

— 

P0 4 

Coliform 
10 9 MPN/ac/day 

Low Res 

.12 

.09 

.002 

.0002 

.0004 

.60 

LM Res 

.43 

.16 

.004 

.0008 

.0006 

.62 

Med Res 

.45 

.18 

.004 

.0008 

.0006 

.63 

High Res 

3.10 

.99 

.006 

.0006 

.0020 

4.9 

Coml 

3.60 

1.35 

.022 

.0060 

.0040 

4.5 

Ind 

6.00 

2.25 

.020 

.0055 

.0030 

5.0 

Agr 

7.20 

2.70 

.001 

.0012 

.00002 

.25 

Pasture 

.24 

.10 

.001 

.0002 

.0002 

.50 

HI Coml 

3.90 

1.44 

.016 

.0065 

.0048 

5.0 

Instl 

3.10 

1.17 

.006 

.0006 

.0020 

6.0 

Natural 

.10 

.04 

.001 

.0001 

.000002 

.0005 

Open 

.24 

.08 

.001 

.0002 

.0002 

.50 
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Dry weather sewage flow was simulated for those basins with significant 
urban land use. Dry weather flow option three was used since it allows 
computations to be made on the basis of land use and population. The 
coefficients used are shown in Table V-2. Domestic, Commercial and 
Industrial coefficients were taken from References 8 and 9, with some 
minor modifications. Pipe Infiltration coefficients were estimated so 
that the quality concentrations did not exceed those for baseflow from 
non-urban subbasins. The coefficients were assumed to remain constant 
for the alternative future. 


TABLE V-2 

Dry-Weather Flow Coefficients 



Domestic 

Commercial 

Industrial 

Infiltration 

low 

|mgd/acre) 

lOO^ 

.03 

.01 

.0005 

Suspended Solids 
(lb/day/ac) 

1.3^ 

1.9 

2.6 

.08 

iettleable Solids 
(lb/day/ac) 

.22- / 

.33 

.44 

.008 

10° 5 

(lb/day/ac) 

.20^ 

.30 

.40 

.002 

1 

(lb/day/ac) 

.04—^ 

.05 

.06 

.0012 

*4 

(lb/day/ac) 

.oi y 

.012 

.02 

.0004 

Collform 

(lO^MPN/day/ac) 

.64^/ 

.86 

.86 

.0125 


1/ gallons/day/caplta for Domestic 
2J pounds/day/capita for Domestic 
3/ 109 MPN/day/capIta for Domestic 


















Since the dry weather flow algorithm In STORM calculates loads and con¬ 
centrations of raw waste water, reductions must be made to account for 
treatment that exists In the study area. An assumption was made that the 
level of treatment remains constant for the alternative future. The 
following removal efficiencies were used for each subbasin having dry weather 
flow. 


Treatment Efficiencies Used in STORM 
(percent) 


Suspended Solids 

87 

bod 5 

87 

Orthophosphate 

80 

Settleable Solids 

87 

Nitrogen 

80 

Collform 

92 


Table V-3 shows predicted tributary loads in pounds for the 10 month 
period of January 1970 through October 1970. While these loads cannot 
be used as an evaluation objective In themselves, they are useful to 
distinguish trends. In every case the predicted loads for 1990C land use 
pattern exceeded those for existing conditions. These loads were not used 
for the Instream analysis. The receiving water analysis was accomplished 
using hourly loads and concentrations for the month of October 1970. 

Storm water guantity and guality were also simulated for the Pendergrass 
detailed study area, however since a receiving water analysis was not 
performed there was no need to predict individual subbasin loadings. 

Table V-4 summarizes the predicted storm runoff guality loadings for the 
Pendergrass detailed study area. 
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TABLE V-3 

Predicted Washoff Of Pollutants 
January through October 1970 
For Athens, Georgia 


Sub- Land Use 
Basin Condition 


Suspended Settleable BODc N PO* Coll form 

(lbs) (lbs) (lbs) (lbs) (lbs) (109 MPN) 

























TABLE V-3 (Cont) 


Predicted Washoff Of Pollutants 
January through October 1970 
For Athens, 6eorq1a 


Sub 

Basin 

Land Use 
Condition 

Suspended Settleable 
(lbs) (lbs) 

BOD, 

(lbf) 

N 

(lbs) 

P0 4 

(lbs) 

Collform 
(109 MPN) 

10 

Exist 

289,780 

37,070 

17,960 

4,613 

1,344 

172,652 


1990C 

333,103 

42,923 

20,516 

5,268 

1,565 

246,885 

11 

Exist 

32,252 

3,586 

1,516 

364 

102 

30,428 


1990C 

41,906 

5,182 

5,482 

604 

185 

84,808 

12A 

Exist 

117,069 

14,402 

7,755 

1,898 

605 

226,020 


1990C 

207,271 

25,616 

13,865 

3,401 

1,161 

504,939 

12B 

Exist 

751,269 

145,068 

49,750 

12,276 

3,920 

1,089,028 


1990C 

735,841 

144,459 

48,679 

11,989 

3,857 

1,144,164 

12C 

Exist 

229,663 

30,845 

17,105 

4,239 

1,216 

439,079 


1990C 

686,401 

102,465 

48,927 

12,159 

4,272 

1,001,348 

13 

Exist 

106,697 

14,890 

9,531 

2,934 

445 

175,192 


1990C 

200,750 

27,880 

19,367 

5,209 

1,792 

458,203 

14 

Exist 

237,387 

29,749 

14,208 

3,649 

1,060 

128,786 


1990C 

251,051 

31,778 

15,142 

3,893 

1,143 

149,505 

15 

Exist 

16,468 

1,982 

873 

197 

52 

14,138 


1990C 

23,865 

2,950 

1,291 

297 

85 

27,180 

16 

Exist 

235,374 

30,904 

13,910 

4,540 

928 

171,134 


1990C 

508,469 

64,824 

33,360 

9,468 

2,571 

872,182 










TABLE V-4 


Storm Runoff Quality Loadings 
Pendergrass Study Area Jan - Oct 1970 



Suspended 
Solids (lb) 

Settleable 
Solids(lb] 

b°d 5 

Ob? 

Nitrogen 

Ob) 

PO 4 

Ob) 

Col 1 form 
(10 9 MPN) 

Existing 

195,244 

23,730 

12,808 

3,099 

933 

214,766 

1990 B 

1,113,700 

142,800 

74,755 

18,467 

6,234 

2,139,334 
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LAND SURFACE EROSION 


Land surface erosion yield was computed by the Universal Soil Loss Equation. 

The equation, as implemented in STORM, requires a continuous hourly precipitation 
record to serve as the prime mover In the analysis. The period of January 
1948-December 1972 at the Athens gage for the Athens area (Winder gage for 
Pendergrass) was used for both existing and future conditions. The K, LS, C, 

P and SDR terms in the equation are shown in Table V-5 [10, 11]. 

The average annual land surface erosion was computed for the period of record 
and from several trial runs It was determined that, for the period of January 
1, 1961 to December 31, 1961 the land surface erosion approximated the average 
annual land surface erosion for the period of record. This shorter period 
was used in subsequent runs to calculate the average annual land surface 
erosion for the various subbasins in the Athens Study Area. 

An Important consideration in the land surface erosion analysis was the 
effect of exposed soil in areas under development. For each grid cell the 
land use for existing and future conditions were compared and the number of 
cells with changed land use were counted. It was then assumed that the 
change In land use will be uniformly distributed over the 15 year period 
1976-1990. Therefore, the area under development for any one year is 
approximately 7% of the total change during the 15 years. For that area 
under development the factors representing the soil cover were modified 
to reflect denuded soil. Specifically, the Cover Factor and Erosion 
Control Factor were set to 100. 
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TABLE V-5 


Soil Erodibility Factors 


SOIL 

NO. 

SOIL 

CODE 

SERIES 

NAME 

K 

SOIL 

NO. 

■a 

SERIES 

NAME 

a 

2 

Ak 



39 

Ln 

Loulsburg 

0.24 

3 

Am 



41 

Me 



4 

An 

Appling 

0.32 

43 

Mg 

Madison 

0.32 

5 

As 



44 

Mi 

Madison 

0.32 

6 

Ax 

Appling 

0.32 

45 

Mm 

Madison 

0.32 

7 

Bfs 

Buncombe 

0.17 


Mm 

Louisa 

0.28 

8 

Ca 



47 

Mv 

Musella 

0.28 

9 

Cb 

Cecil 

0.32 

48 

Pa 



10 

Ce 



49 

Pf 

Pacolet 

0.32 

11 

Cf 



50 

Pg 

Pacolet 

0.32 

12 

Ci 

Colfax 


51 

Ph 

Pacolet 

0.32 

13 

Coa 

Congaree 


52 

Pi 

Pacolet 

0.32 

14 

Cob 

Chewacla 


54 

Pt 



17 

Cy 



55 

Rc 



18 

CY 

Cecil 

0.32 

56 

Rok 

Rock 

0.00 

19 

CZ 

Cecil 

0.32 

57 

Tf 



21 

Oh 

Davidson 

0.32 

58 

To 



22 

Dq 

Davidson 

0.32 

61 

Wq 



23 

EW 



62 

Wk 

Worsham 


25 

Ge 



63 

Wos 

Wehadkee 


29 

Gr 



65 

LD 

Louisburq 

0.25 

31 

He 







33 

H1 1 








NOTE: SOILS WITH NO "K” VALUE USED THE DEFAULT OF 0.32 
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TABLE V-5 


< 


Length for LS Factor 


SLOPE * 


LENGTH (ft) 


0 - 2.00 200 

2.01 - 6.00 275 

6.01 - 10.00 175 

10.01 - 15.00 75 

15.01 - 25.00 50 


Cover and Erosion Control Factors 


Land 

Use 

Cover Factor (C) 

Erosion Control 
Factor % (P) 

1 . 

Open 

1.3 

05 

2 . 

Low Residential 

0.3 

85 

3. 

Medium Residential 

0.3 

70 

4. 

High Residential 

0.3 

60 

c 

J • 

Agricultural 

40.0 

95 

6 . 

Industrial 

10.0 

40 

7. 

Commercial 

1.2 

20 

8 . 

Pasture 

2.0 

95 

9. 

High Commercial 

1.0 

10 

10 . 

Institutional 

10.0 

40 

11 . 

Roads 

5.0 

60 

12 . 

Natural 

0.3 

95 









TABLE V-5 


Sediment Delivery Ratio (SDR) 


WATERSHED NO.* 

AREA, ac 

SDR 

1 

6669 

0.18 

2 

10444 

0.16 

3 

2272 

0.23 

4 

7915 

0.17 

5 

7648 

0.17 

6 

16301 

0.14 

7 

6101 

0.18 

8 

12164 

0.16 

9 

6242 

0.18 

10 

1290 

0.25 

11 

5946 

0.19 

12 

3460 

0.21 

13 

3438 

0.21 

14 

5039 

0.19 

15 

3103 

0.21 

16 

11253 

0.16 

Pendergrass 

7067 

0.17 











All other factors In the soil loss equation remained the same as In 
the developed condition. The predicted land surface erosion for the 
Athens and Pendergrass study areas are shown In Tables V-6 and V-7. 
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Table V-6 


Athens* Georgia 


Average Annual 

Land Surface Erosion (tons) 

Watershed No 

.* Existing 

1990 C 

1 

76900 

72000 

2 

112000 

115300 

3 

5000 

6300 

4 

85000 

83900 

5 

51100 

55700 

6 

118700 

306200 

7 

10100 

16900 

8 

222700 

219900 

9 

93100 

88500 

10 

17900 

17100 

11 

96800 

93800 

12 

24500 

24400 

13 

23500 

23100 

14 

89000 

67000 

15 

4800 

5200 

16 

163100 

157400 


♦Savannah District Watershed Identification 
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TABLE V-7 


Average Annual 
Land Surface Erosion 
Pendergrass Area 


Existing Land Use 105220 tons 

Alternative B Land Use 94910 tons 
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RECEIVING WATER 


Analysis of Existing Condition 

The WQRRS model accepts Input tributary conditions derived using the STORM 
model on each of the 32 subbasins (see Figures III-2, III-3 and III-4) 
draining Into the portion of the Oconee River within the specified study 
boundaries and Imposes these loadings on a base flow condition. The two 
sewage treatment plants and the Athens water Intake are accounted for based 
on mean monthly data from the State of Georgia, except for unmeasured 
parameters which were then estimated from textbook average conditions [12]. 
These Input data are shown In Table V-8. An accounting Is made of the mass 
balance at each tributary location and the resulting mixture Is transferred 
(l.e., routed) downstream to the next tributary location with the proper 
reactions and Interactions being calculated according to the estimated 
travel time between tributaries and the Input system coefficients shown 
In Table V-9. 

The Initial quality condition for selected locations Is shown In Table V-10. 
The quality at all other locations Is obtained by linear Interpolation. The 
values shown for river mile 34.8 and 31.5 on the North Oconee and Middle 
Oconee respectively are also the base flow quality conditions which enter 
the study area at the upper boundaries. 

The base flow condition on the tributaries during non-storm periods Is 
dependent on the proportion of the drainage area having residental land 
use. Table V-U shows the tributary base flow used as Inflow to the Oconee 
during non-storm periods. 
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TABLE V-8 

INPUT DATA FOR SEWAGE TREATMENT PLANTS AND 
WATER TREATMENT PLANTS FOR EXISTING AND 
ALTERNATIVE FUTURE C LAND USE 


Parameter 1/ 

North Oconee 
STP(R.M. 19.9) 

Exist " Alt.C 

Middle Oconee 
STPfR.M. 17.3) 
Exist JCTt.C 

Athens Water 
Intake(R.M. 23.7) 
Exist Alt.C 

Q (cfs) y 

10.1 

15.3 

3.25 

5.0 

16.6 

24.9 

Temperature (°C) 

3 / 

3/ 

3/ 

If 



DO (assumed) 

0 

0 

0 

0 



BODS 4/ 

99 

99 

84 

84 



Collform (assumed) 
(MPN/100 ml) 

ZOO 

200 

200 

200 



Detritus 

(25% of susp. solids) 

8.25 

8.25 

7.25 

7.25 



NH3 [5] 

10 

10 

10 

10 



N03 [5] 

20 

20 

20 

20 



N02 [5] 

.05 

.05 

.05 

.05 



P04 [5] 

12 

12 

12 

12 



TDS 

244 

244 

154 

154 



Algae (assumed) 

.001 

.001 

.001 

.001 



Zooplankton 

(assumed) 

.001 

.001 

.001 

.001 



pH (units) 

6.7 

6.7 

7.3 

7.3 



Akallnlty (assumed) 

100 

100 

100 

100 




1/ mg/1 except as noted. 

y Flow for Alternative C equals existing flow times estimated.proportional Increase 
In population. 

y Water temperature equals mean dally air temperature minus 2*C [13] except during 
storm events when water temperature equals the hourly air temperature. 

4/ Uncorrected for NHj, N0 2 and Detritus oxygen demand Cl 3 3 • 
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TABLE V-9 

INPUT SYSTEM COEFFICIENTS 

REACTION PATF MULTIPLIER PARAMETERS 


CALIBRATION MAGNITUDES CALIBRATION TEMPERATURES 



Kl 

K2 

K3 

K4 

T 1 

T 2 

T 3 

T 4 

ALGAE 1 

. 10 

.98 

.98 

. 10 

5.0 

22.0 

25.0 

34.0 

ALGAE 2 

.10 

.98 

« 9 R 

.10 

10.0 

28.0 

30.0 

40.0 

ZOOPLANKTON 

.10 

.98 

.98 

.10 

5.0 

28.0 

30.0 

39.0 

BENTHIC ANIMALS 

. 10 

.9e 

.98 

. 10 

5.0 

22.0 

25.0 

33.0 

FISH 1 

.10 

.98 

.98 

. 10 

5.0 

20.0 

20.0 

25.0 

FISH 2 

.10 

.9 8 

.98 

.10 

10.0 

27.0 

30.0 

33.0 

FISH 3 

. 10 

.98 

.98 

.10 

5.0 

22.0 

30.0 

36.0 

BOD 

.10 

.98 



4.0 

30.0 



NH 3-N 

.10 

.98 



4.0 

30.0 



N02-N 

.10 

.98 



4.0 

30.0 



DETRITUS 

.10 

.98 



4.0 

30.0 




DECAY CCJEFf ICIENTS, PER DAY 

MAX VALUE 

BOD 

. 100 

NH3-N 

.050 

N02-N 

. 2C0 

DETRITUS 

.001 

COL I FOP M (AT 20 OEG C) 

.500 


010 TEMPERATURE COEFFICIENT FCR CCLIFCRM 1.040 

CHEMICAL COMPOSITIONS OF BIOTA 

C N P 


ALGAE 

. 500 

.090 

.012 

ZOOPLANKTON 

.500 

.090 

.012 

FISH 

.500 

.090 

.012 

BENTHOS 

.500 

.090 

.012 

detritus 

.500 

.090 

.012 



TABLE V-9 (cont'd) 


RESPIRATION PATES,PFP DAY 

MAX VALUE 




PHYTOPLANKTON 

.500E-01 




ZOOPLANKTON 

.200E-01 




FISH 

.I00E-02 




BENTHOS 

.10CF-02 




DETRITUS SETTLING, METERS/OAY 

.15000 




OTHER PHYTOPLANKTON DATA 





SETTLING, METER/DAY 

.15000 

.15000 



OXYGENATION FACTOR 

1.600 




PRFFERANCE 

.670 

.330 



SELFSHADING PEP “G/L/M 

0 




MAXIMUM SPECIFIC GROWTH PATE, 

PER DAY 




PHYTOPLANKTON, 2 GROUPS 

. 1 COF 01 

• 200F+-01 



ZOOPLANKTON 

.15CE+00 




FISH, 3 GROUPS 

•200E-01 

.250E-01 

. 200E-01 


BENTHOS 

.2C0F-01 




HALF-SATURATION CONSTANTS CF 

ALGAC 





LIGHT 

C n 2 

N 

PJ4 

ALGAE L 

.003 

. 020 

.200 

. 030 

ALGAE 2 

.005 

.020 

.100 

.050 


HALF-SATURATIGN CONSTANTS for 

ZOO, FISH AND BENT HO 

ZOO GRAZE ON ALGAF 

.550 

FISH 1 GPAZE GN ZOO 

. 050 

FISH 2 GRAZE ON ZOO 

.050 

FISH GRAZE ON RENTmCS 

500.000 

BFNTHOS GRAZE ON SFD*T 

50.C 00 


STOICHIOMETRIC EQUIVALENCE 

02/NH3 
C2/N02 
02/0E TRITUS 
02/BIOMASS 
C02/B00 


OF CHfviCAL TRANSFORMATION 

3.300 
I .200 
2.000 
2.000 
.200 
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TABLE V-10 

INITIAL QUALITY CONDITION 


Location 

Parameter 

Magnitude 1/ 

Source 

North Oconee-RM 34.8 

BODS 

.5 

Smith [5] 


Detritus 

5 



Sediment 

(gm/m2) 

16 



O 

Benthos (gm/m ) 

.9 



NH- 

no 3 

.03 

.24 

State of Georgia gage 
at Athens Intake, 
average 1974 


NOp 

.01 



P0 4 

.03 



pH (pH units) 

7.1 



Alkalinity 

36 



Coliform 

430 

Athens Intake 


(MPN/100 ml) 


10/24/74 


Temp. (°C) 

f(alr)—^ 

Willey A Huff Tl3] 


DO 

8 

Assumed at 80^ of 




D0 sat at 15°C 


TDS 

100 

Assumed 


Algae 

.001 



Zooplankton 

.001 



Fish 1 (Kg/mi) 

10 

Smith [5] 


Fish 2 (Kg/ml) 

30 



Fish 3 (Kg/ml) 

40 



60 













TABLE V-10(cont'd) 
INITIAL QUALITY CONDITION 





>^rv. 


OP fc**- 














TABLE V-10 (cont'd) 
INITIAL QUALITY CONDITION 


Location 

Parameter 

Magnitude 1/ 

Source 

Middle Oconee-RM 31.5 

B0D5 

.5 

Smith [5] 


Detritus 

5 



Sediment 

(gm/m 2 ) 

16 



p 

Benthos (gm/m ) 

.9 



NH, 

.10 

USGS gage on Middle 



.60 

Oconee NH 3 : 2/09/70 


1U 3 

NO 3 + NO 2 : avg of 


NO, 

.01 

11/18/70 and 5/26/70 


PO 

.06 

PO 4 : 9/02/70 pH and 


KU 4 

Coliform: 11/18/70 


pH (pH units) 

7.5 



Coliform 
(MPN/100 ml) 

930 



Alkalinity 

25 

Alkalinity: avg of 
9/02/70 and 11/18/70 


Temp. (°C) 

f(air)^ 



DO 

8 

80% of D0 sat at 15°C 


TDS 

100 

Assumed 


Algae 

.001 



Zooplankton 

.001 



Fish 1 (kg/ml) 

10 

Smith [5] 


Fish 2 (kg/mi) 

30 



Fish 3 (kg/mi) 

40 



6 2 















TABLE V-10 (cont’d) 
INITIAL QUALITY CONDITION 


Location 

Parameter 

Magnitude 1 J 

Source 

Middle Oconee-RM 17.5 

B0D5 

NH- 

»°3 

*>4 

Collform 

All other parame 

.5 

.08 

.22 

.10 

660 

:ers same as R.M. 31. 

STORM Base Flow 

5 

Middle Oconee-RM 17.0 

Detritus 

(gm/m 2 ) 

Sediment 

(gm/m 2 ) 

All other parame 

5.5 

18.4 

;ers same as R.M. 17. 

Smith [5] 

5 

Main Oconee-RM 0.0 

All parameters a 

*e same as Middle Occ 

nee R.M. 17.0 


1 / mg/1 except as noted 

2 / water temperature equals hourly air temperature during a storm and equals mean 
daily air temperature minus 2°C during non-storm periods [13]. 
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The Initial and base flow quality conditions are arbitrarily accepted base 
conditions since essentially no gaged data exists for the study period. 

All the final results must be Interpreted relative to this base condition 
since most of the water quality calculations are non-linear (l.e., effect of 
saturation values and temperature corrections on all reaction rates). 

Table V-12 shows an example of a statistical summary of the water quality 
condition at a random point along the river for gristing land use. The 
critical values (l.e., maximum or minimum) for some of these parameters at 
various locations along the river have been plotted in Figures V-l to V-6 
to show a river water quality profile for the most critical condition 
occurring during the period simulated for dissolved oxygen (DO), ammonia 
(NHj), nitrate (N0 3 ), phosphate (PO^), log coliform bacteria and 5-day 
carbonaceous biochemical oxyqen demand (B0D 5 ). 

The existing water quality condition seems to meet all the Georgia State 
Water Quality Standards (i.e.. Table V-l3) except for coliform bacteria 
which may exceed the standards 5-10% of the time, from river mile 25 to 
12.3 respectively on the North Oconee and throughout the study length of 
the Middle and Main Oconee about 10-16% of the time. 

The impact on the water quality in the Oconee due to the various tributaries 
and sewage treatment plants is shown in Figures V-l to V-6. Major point 
source Impacts on the North Oconee River are summarized in Table V-14, and 
on the Middle and Main Oconee River in Table V-l5. 
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Remarks in Tables V-14 and V-15 concerning nutrients having significant 
Impact refer to the potential impact on algae production in non-turbid 
water. Unless significant improvement occurs in the turbidity of the 
Oconee River, this potential will be not be realized. 

Sample graphical results of the simulations are shown in Appendix B. 


i 











TABLE V-12 

WATER QUALITY AT RIVER MILE 11.5 
EXISTING UNO USE 


POST-PHOCESSuW FUR -awRS APRIL 1976 

MYOROLOGIC EmGINFFHIMG CENTER DAVIS, CA 


ucgnee river water 'Jijalitv study *»wqr«s statistical post-processor** 

REACH a MIOOLE UCUNEf RIVER (R.m. 17,0-8.01_ 

QUALITY data RaSLC' Ur. ExISTTnG UAiwu USE 

• ***•«*•*•*•**«*••• I»,PUT DATA MtMtMUAMMtttMM 


BEGINNING OF REACH RtvEL n tLE 17,00 

END OF REACH RIVER MILE _8,00 _ 

SURREACH LENGTH (MILES) ,50 

COMPUTATION INTERVAL (HOURS) _ - 2- _ _ 

-A-JBST- OA-Y- C1F —StMuIA. C m» L-gE JUi lQ-274 t t OCT 701— 

LAST day of simulation RERinO 304 (31 uCT 70) 

NUMBER OF days in simulation period 31 . ... . 


OBSERVATIONS at RI 
first day flF STunr 
-LAS T. PAY- -OF SLUOY^ 
NUMBER OF DAYS In 

v£R mu f 

PFRInn 

. - . 

tl.50 

274 ( l 

_30-4_Ui 

UCT 70) 
lict /m 

-*-- 

study otMlf'O 


31 


on. nF 

rater duality parameters at river mile 
number OF SIMULATION points 

11.50 

373 

E A X A A AAA — . - _ 



- simulation 

VALUES 


(SIMULATED-OBS,) 

OBSERVED 

parameter 

minimum 

“a Xjmum 

MEAN 

SID.DEV. 

mean std.dev. 

. VALUES 

FLOw 

8.9 

37.1 

12.7 

5.9 



TEMP 

9.2 

24.9 

17.7 

3,4 

0,0 0,0 

0 

__AX Y 

_4.0 

U.2 . . 

_9.5- 

.6- 

0,0 0,0 

0 

NH3 

,0 80 

.475 

,3Rl 

.084 

0.000 0,000 

0 

Nl) J 

.220 

1,215 

1.063 

.166 

O.OOU 0,000 

0 

POa 

.100 

.566 

,440 

.103 

0,000 0.000 

0 

ALKA 

2b.o 

3) .2 

30.0 

.6 

0.0 0,0 

0 

LOG COL I 

2.78 

4.44 

3.03 

.“2 

0.00 0.00 

0 

—TOS - 

-LQO-, 

Lai.- 

.—1 r — 


0, 0, 

A 

PH 

7.3 

7.7 

7.6 

. i 

o.o 0,6 

0 

ROD 

• 5 

n.t 

2.9 

.7 

0.0 _ 0,0 

0 



Parameter 


*ATER UOalITy POINTS EXCEEDING STaNPARD 

_STAIDARD . _-NUMBER- - - PERCENT._ 


TEMP 

OXY 

LOG CULX 
PM 
PH 


32,2 MAX. 

4,0 MIN, 
3,6'* “AX, 
6 ,'l “IN. 
-8,5 MAX. 


0 

0,00 

0 

0.00 

52 

13.94 

0 

0.00 

0 

0.00 



















TABLE V-13 

GEORGIA WATER QUALITY STANDARDS 
AND 

OCONEE RIVER STREAM CLASSIFICATION 


DRINKING WATER STANDARDS 
(waters requiring treatment) 


Coliform 

maximum 

<«nQn MPN/inn nl 

Dissolved Oxygen 

minimum 

4 mg/1 

(warm water fish) 


pH 

minimum 

6.0 


maximum 

A. 5 

Temperature 

maximum 

90°F 

NO^ as Nitrogen 

maximum 

10 mg/1 


FISHING WATER STANDARDS 


Coliform 

maximum 

4QOH MPN/lon ml 

Dissolved Oxygen 

minimum 

4 ma/1 

(warm water fish) 


pH 

minimum 

6.0 


maximum 

8.5 

Temperature 

maximum 

00° F 


STREAM CLASSIFICATION 



North Oconee 

River 

R.M. 

R.M. 

34.8-23.7 

23.7-12.3 

Drinking Water 
Fishing Water 

Middle Oconee 

River 

R.M. 

31.5-21.8 

Drinking Water 

Middle A Main 

Oconee River 

R.M. 

21 .8- 0.0 

Fishing Water 








NORTH OCONEE 



RIVER MILE 





NORTH OCONEE 



Figure V-2. Existing Land Use 
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RIVER MILE 
















Figure V-4. Existing Land Use 
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Figure V-5. Existing Lond Use 
















NORTH OCONEE 



Figure V-6. Existing Land Use 





















TABLE V-14 

POINT SOURCE IMPACTS ON THE NORTH OCONEE RIVER 
EXISTING LAND USE 


RIVER 

MILE 

PROBABLE 

POLLUTANT 

SOURCE 

PARAMETER 

IMPACT 

-. 

REMARKS 

26.0-23.7 

Unknown 

DO 

•B mg /1 

minor 


Subbasin 1A * 2 

Coliform 

9000 MPN/100 ml 

significant 

23.2-21.7 

Subbasin 3 

bod 5 

2.5 mg /1 

mi nor 



nh 3 

.2 mg /1 

significant 



no 3 

.5 

significant 



PD 4 

.2 mg /1 

significant 

19.9-19.3 

STP 

DO 

.6 mq /1 

minor 



BOD 5 

4 mg/1 

minor 



nh 3 

.6 mg /1 

sianificant 



N 0 3 

1.1 mg /1 

significant 



P0 4 

.7 mg/1 

significant 


TABLE V-15 

POINT SOURCE IMPACTS ON THE MIODLE AND MAIN OCONEE RIVER 

EXISTING LAND USE 


m 

PROBABLE 

POLLUTANT 

SOURCE 

PARAMETER 

IMPACT 

REMARKS 

22.5-17.0 

STP * Sub- 

NH, 

.1 mq/1 

significant 


basin 6A R 6B 

no 3 

.2 mq/1 

significant 



P0 4 

.1 mn/1 

significant 



BOD 5 

1.6 mq/1 

minor 

12.3 

North Oconee 

NH 

.1 mg/1 

significant 



"4 

.1 mg/1 

significant 

7.2 

Subbasin 16 

BOD 5 

.5 mo /1 

minor 



nh 3 

.1 mg /1 

significant 
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Analysis of Alternative Future C 


The STORM results derived from alternative future C land use condition 
were imposed on the same initial river quality condition as defined for 
existing land use in Table V-10. The tributary base flow quality condition 
used for alternative C land use is shown in Table V-16. Table V-17 shows 
an example of a statistical summary of the water quality condition at a 
random point along the river for alternative future C land use condition. 
The critical values for some of these parameters at various locations 
along the river have been plotted in Figures V-7 to V-12 to show a river 
water quality profile for the most critical condition occurring during 
the period simulated for dissolved oxygen (DO), ammonia (NH 3 ), nitrate 
(NO-j), phosphate (PO^), log coliform bacteria and 5-day carbonaceous 
biochemical oxygen demand (BOD^). 

The water quality condition for alternative future C land use condition 
exceeds the Georgia State Water Quality Standards (i.e.. Table V-13) 
similar to the existing condition. 

Sample graphical results of the simulations are shown in Appendix R. 
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TABLE V-16 

TRIBUTARY BASE FLOW QUALITY DATA]/ 
ALTERNATIVE C LAND USE 
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TABLE ¥-17 

WATER QUALITY AT RIVER MILE 11.5 
ALTERNATIVE FUTURE C LAND USE 


*»* ••«**«*«« *»» * ******•«.«»*»,*,** A *** AA____ 

POST-PRi’CFSSLR FIR -iQRRS APRIL 1576 
H*D«OLOGtC ENGINEERING CFNTER DAVIS. C A _ 


OCl>EF RIvER -ATER Quality study aahnRHS STATISTICAL P03T-PR0CES30H*» 

HEACf 4 a "lO'JLt. QCUREE RlV£R (rt.M. 17*0-6.0) _ _ 

QUALITY DATA RASED Of; ALTERNATIVE C LAND USE 


INPUT QATA 


BEGlNNtuG Of reach RIVER RILE 
END OE REACH RIVER H ILE 

SUflHEACM LENGTH (M t lE S) 

COMPUTATION TNTtRVAL (HOURS) 

FIRST-QA.Y- JOE. SlHOl_AXLUil_P.tilICC_ 2J.U t l__UC.I— Z&l 

L*ST Oav (jF 31 **-JL ATI "IN pfKluO 304 (31 UC T 70) 

NUHhEW nF DAYS 1*1 SlMjLATIf’M PERIOD 31 

OBsEnVATti.lNS AT RIVER “H£ 11.so 

FIRST Day of STUn> PERIOD 27“ ( 1 UCT 70) 

_lasi daj__£le_ iLuaY -PiauDa___sou u i _uc Lmi 

NU”flER llF DAYS In STUDY PfcRIi'O 31 


17.00 

a .oo 

.50 

2 


•iateR quality parameters at river mile 

NUMBER C'F SIMULATION POINTS 

11.50 

373 

—EHBflfi 

OBS.) 

-.DEV. 

0.0 
n, o 

-NO»_Qf _ 

urserved 

values 

0 

(L 

parameter 

FLO* 

V E M P 

OXY 

I- INI HU- 

5.7 
q .2 
_ _a.u 

- SI -Ul A T JPN 
MAXIMUM. 

37.4 

2o. 6 

11.2 

VALUES --- 

H£AN -STD,DEV, 

13. 4 6.2 

17.6 3.4 

(3IMULATE0- 
-..MfcAN STo 

0,0 

NH 3 

. Ofl i) 

.625 

.51 *) 

.115 

0.000 

0,000 

0 

NO 3 

.221' 

t .531 

1.326 

.230 

o.ooa 

0.000 

0 

PUU 

. 10(1 

.740 

.5^0 

.135 

0.000 

0,000 

0 

4LKA 

25..) 

32. I 

3o.« 

1.0 

0,0 

0.0 

0 

LOG CUl.I 

2 . *2 

u.53 

3.10 

.“6 

0,00 

0.00 

0 

T03 - 

_loa. 

— lat.— 

tao. — 

I. - 

-- 

-a. 

<L 

PH 

7.3 

7.6 

7.5 

. 1 
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WOO 

.5 

5.1 

3.5 

.5 

0,0 

0,0 

0 



-ATE'S QUALITY 

POINTS EXCEEDING STANDARD 


PARiMETt R 

STANDARD 

number _ 

-PERCENT 
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TEMP 

32,2 MAT, 
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0,00 


oxy 

4,0 HIN. 
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NORTH OCONEE 



Figure V-7. Comparison of Water Quality Due to Changing Land Use 
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V-8. Comparison of Water Quality Oue to Changing Land Use 
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Figure V-9. Comparison of Water Quality Due to Changing Land Use 
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Figure V-10. Comparison of Water Quality Due to Changing Land Use 















Comparison of Water Quality Due to Changing Land Use 
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-12. Comparison of Water Quality Due to Changing Land Use 




















Impact of Alternative Future C 

The water quality Impact of alternative future C land use is shown in 
Figures V-7 to V-12. They have been summarized in Table V-18. 


TABLE V-18 

WATER QUALITY IMPACTS DUE TO CHANGING FROM 
EXISTING LAND USE TO ALTERNATIVE FUTURE C 



" MaGnItud 
( mg/1) 

E LOCATION 

(river mile) 

S 


BOD, 

1-2 

North Oconee 

25-12.3 

minor 

D 

1 

Middle Oconee 

25-0 

minor 

DO 

.2-.4 

North Oconee 

22-12.3 

minor 


.2 

Middle Oconee 

18-15 

minor 

NH, 

. 1-.2 

North Oconee 

25-12.3 

significant 


.1-.2 

Middle Oconee 

25-0 

significant 

NO, 

. 1- .4 

North Oconee 

25-12.3 

significant 

0 

.1-.3 

Middle Oconee 

25-0 

significant 

PO, 

.1-.3 

North Oconee 

25-12.3 

sionificant 


.1- .2 

Middle Oconee 

25-0 

significant 

Coliform 

(MPN/lOOmT) 

3,000- 

10,000 

North Oconee 

33-12.3 

significant 


3,000- 

10,000 

Middle Oconee 

25-n 

significant 


Remarks in Table V-18 concerning nutrients having sinnificant impacts refer 
to the potential impact on algal production in non-turbid water. Unless 
significant improvement occurs in the turbidity of the Oconee River, this 
potential will not be realized. 

In general, the sources of the increased pollutants due to changing land use 
are the same as those defined f or existing conditions in Tables V-14 and V-15. 
Concentrations of pesticides, heaw metals and other parameters not speci¬ 
fically mentioned were not evaluated in this study. 









GRID CELL SEDIMENT TRANSPORT INVESTIGATIONS 


Introduction 

A distributed parameter, structure imitating model was developed for the 
calculation of land surface erosion and deposition. The phenomena 
simulated in the model are: rainfall-runoff, runoff accumulation and 
distribution, detachment of soil by rainfall, transport of detached soil 
by runoff, scour by runoff, and deposition. Application of the model to 
laboratory test data yielded encouraging results. The application to a 
watershed in this study was unsuccessful, however, because of unsatis¬ 
factory topographical information. 

Model Description 

The model performs calculations on a cell-by-cell basis. The direction 
and velocity of runoff are determined from topographical information 
Imbedded In the data base. A steady state process Is assumed. Details 
of the computations are given below. 

Rainfall-Runoff 

The very simple "rational formula" was used to generate runoff from each 
cell. The runoff from any cell Is: 

Q * ciA 
where: 

0 * discharge, cfs 
c * runoff coefficient 
1 * rainfall intensity, in/hr 
A * cell area, acres 
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The value of the coefficient, c, was determined from a combination or 
hydrologic soil type and land use as shown in Table V-19. The values 
in the table are judgemental and have not been calibrated. 

It is believed that use of this rainfall-runoff relationship is justified 
because of the small scale of the cells and the steady state nature of 
the simulation. The method for accumulating runoff from individual cells 
Is described below. 

Runoff Accumulation 

Consider a typical cell (I, J) and its eight neighboring cells as shown 
in Fig. V-13. If any of the neighboring cells are at higher elevations, 
a portion of runoff generated at those cells will reach cell (I, J). The 
runoff generated within cell (I, J) is found by the rational formula and 
added to the sum of all the contributions from higher cells to give the 
total discharge passing out of cell (I, J). This discharge Is evenly 
distributed among all neighboring cells of lower elevation. For this 
reason computations must proceed from higher elevations to lower which 
requires that the data bank first be sorted by elevation. Note that, 
since the sediment moves with the runoff, this portion of the calculations 
also determines the paths that the sediment takes. No runoff (or sediment) 

Is passed between cells of equal elevation. 

Soil Detachment by Rainfall 

Following a suggestion by Foster and Meyer [14] it is assumed that rate of 
detachment Is proportional to the square of the rainfall Intensity with 
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the proportionality constant being the erodlblllty factor (K) In the 
universal soils loss equation. The formula used Is: 

RE = a (K)(Ex)(A)(i 2 ) 
where: 

RE * rate of soil detachment by rainfall, tons/hr 

a * an empirical constant, the value used (0.0002) was based on 
very sparse data 

K * soil erodlblllty factor in the universal soils loss equation, 
related to soil type 

Ex = percent of the cell area exposed to rainfall, related to land 
use as shown on Table V-19 

? 

A = cell area, ft 
i * rainfall intensity, in/hr 

In this study, the rainfall Is assumed to be uniformly distributed over 
the basin and all cells are assumed to be of equal size. Therefore, 
only the values of K and Ex change from cell to cell to reflect spatial 
land use variation. 

Sediment Transport by Runoff 

The hydraulics of the flow must be further defined before sediment transport 
calculations can begin. To define the hydraulics certain Important 
assumptions must be made. The formation of rills and gullies Is Important 
to the runoff hydraulics, however no means of predicting their formation 
or ultimate size was found. Therefore, It was assumed the runoff between 
cells occurs as sheet flow. The width scale, B, of the runoff Is, therefore, 
related to the total length of a cell boundary and the number of neighboring 
cells to which runoff Is passed (those with lower elevations). If flow goes 
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to all eight neighbors, B equals the total length of the cell coundary; 
if flow goes to only one cell, B is one-eighth of the total, etc. From 
the width, B, a typical depth of flow, y for each outflow path is calculated 
from the Manning eq.: 

y , Qn _ 

1.486 B S 1/2 

where: 

y » depth of flow, ft. 

Q * total discharge from cell, cfs 
n = Manning's n 
B = width of outflow path, ft. 

S * slope to particular cell, difference in elevation divided by 
distance between cell centroids 

Sediment transport rates and erosion or deposition by the runoff are 
based on a simple DuBoys relationship: 

TC * CS(t - xm) b 
where: 

TC * transport capacity of any single outflow path, tons/hr 

CS * transport coefficient, related to representative grain size: 

CS 3 52.3 where 0 is grain size in millimeters, 

empirically related to soil type 

2 

t 3 shear stress, Ib/ft 

t * yyS, where y 3 unit weight of water and y and S are as previously 
defined 

2 

rm = critical shear stress below which no transport occurs, Ib/ft 
related to grain size 

b 3 width of an Individual outflow path. Constant at one-eighth 
the total length of a cell boundary 
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Whether scour or deposition occurs alonq any particular outflow path depends 
upon whether the sediment load is less than or greater than the transport 
capacity. The load to any outflow nath is calculated as follows: 
g * (GI + RE) (Q0PS/O0) 
where: 

g * sediment inflow to any particular outflow path, tons/hr 

GI = total rate of sediment inflow to the cell from neiqhborinq cells 
of higher elevation, tons/hr 

RE 3 rate of sediment detachment by raindrop within the cell, ton/hr. 

Q0PS 3 runoff following individual outflow path, cfs, equal to Q0 divided 
by the number of outflow paths 

QP 3 total runoff from cell, cfs 

If g is less than the transport capacity, erosion occurs along the outflow 
path. The actual transport rate for that path is calculated by the 
following: 

G2N 3 q (1 - Ex) + ((1 - F) q + F(TC)Ex) 

where: 

G2N 3 sediment transported out of cell alonq anv given outflow path, 
tons/hr 

F = a relaxation factor if F = 0, outflow = inflow; if F = 1, 
outflow 3 transport capacity 

Other symbols are as previously defined. The exposure factor appears 
because erosion can only occur where the soil is available. 


If g is greater than the transport capacity, deposition occurs and the 
outflowing load is calculated as follows: 









G2N « (1 - F) g + F (TC) 


All symbols have been defined. The exposure does not appear because 
deposition can occur everywhere. 

These values are added to the inflowing load (Gl) of the neighboring 
cells. 
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Application to Laboratory Data 

The algorithm for calculation of land surface erosion and deposition 
described above was tested by comparing calculated erosion rates with 
those measured In a laboratory. While such a test does not constitute 
rigorous verification, it can be used to evaluate the general validity 
of the approach and Identify some inadequacies. 

The laboratory test data used [15] Is from a 5 foot by 16 foot plot. 

It was modeled using 25-1 by 1.6 feet rectangular cells. The rough¬ 
ness value (n a 0.022) and runoff coefficient (C * 0.97) were based on 
measurements made during the experiment. Although several slopes were 
tested, only the 10% slope condition was modeled. This slope was re¬ 
flected In elevations assigned to the Individual cells. The simulated 
rainfall Intensities were used in the program and the calculated weight 
of sediment transported to the bottom of the plot compared with that 
measured. The results are shown in Table V-20. 

The same set of coefficients was used for all rainfall intensities In¬ 
dicating that the functional relationship between rainfall and erosion 
used Is reasonable. A relaxation coefficient, F, of zero had to be used. 

A zero F indicates that the runoff has sufficient transport capacity to 
carry all the sediment produced by raindrop erosion. Apparently the 
transport capacity calculated by the DuBoys relation was too small. 

This could be due to the assumed runoff hydraulics being Inadequate, 
or the critical shear stress not being appropriate for land surface 
erosion. The critical shear stress was taken from a Shield's diagram [16] 
developed for open channel flow. 
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TABLE V-20 


Rainfall 

Intensity 

(In/hrj 

Measured 

Wt. of Sediment 
(lbs/hr) 

Calculated 

Vft. of Sediment 
(Ibs/hr) 

1.25 

5.3 

7.6 

2.25 

27.1 

24 i 

i 

3.65 

67 

66 


4.60 


106 


104 









Application to Sandy Creek Basin 

The model was applied to the lower 16 square miles of the Sandy Creek 
watershed. An existing detailed data bank was available for this area. 

The spatial variation of various parameters within the watershed was 
described using a total of 9208 grid cells. The variables used from the 
data bank were: cell elevation, soil type, hydrologic soil type and land 
use. A topographic map of the basin is shown in Figure V-14. 

Runoff coefficients were determined from land use and hydrologic soil 
type as given in Table V-19. Also shown are the exposure values assigned 
to the various land uses. Descriptions of the land uses are given in 
Chapter III, and soil erodlbility factors (K) in Table V-5. After several 
runs and mapping of computed discharges and sediment loads, a basic data 
problem was identified which prevented completion of the apolicatlon. 

This problem Is discussed in detail below. 

Topographic Data Problems 

The procedure used relies on topographic data in the form of an elevation 
for each cell. Differential elevations between neighboring cells drive the 
runoff calculation. If the cell elevations are truly representative of the 
topography of a drainage basin, every cell but one will have at least one 
outflow path. The exception is the outlet on the watershed boundary. 

The cell elevations In the Oconee study were manually assigned from a base 
topographic map. This procedure resulted in many cells for which no outflow 
path existed. Out of 9208 total cells, 334 had no outflow path. Since no 
runoff or sediment can pass through a cell with no outflow path, the 
calculations for all downhill cells are erroneous. 
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Several "smoothing" algorithms were tried to insure that all cells had 
outflow paths. Of these, the best appeared to be the following: set the 
no outflow cell's elevation equal to the average elevation of the lowest 
and second lowest neighboring cells. This guarantees that the cell in 
question has an outflow path, but may eliminate the only outflow path a 
neighboring cell has. The algorithm worked fairly well; after several 
passes the number of no outflow cells was reduced from 334 to about 18, 
which could not be further reduced by successive application of the 
algorithm. Unfortunately, this was still too high to yield acceptable 
results, as shown on Figure V-15. Happed on that Figure are cell-by-cell 
discharges generated by 1-1nch of rain in one hour. The tendency for the 
runoff to accumulate In stream channels can readily be seen. Note also, 
however, the effect of cells with no outflow path (indicated by large 
dots). The effect Is to disconnect the basin, so that runoff from upper 
portions of the basin does not pass through. Sediment transport calculations 
based on this runoff pattern are, of course, meaningless. 

If a workable algorithm cannot be developed for editing the topographic 
data, the procedure for initially determining cells' elevations must be 
modified. Automatic interpolation is one possibility. This process 
would rarely produce neighboring cells of exactly the same elevation. 

Another procedure being investigated describes topography as an array 
of triangular elements. Elevations are prescribed at the vertices and 
vary linearly within each element. 
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ADAPTION OF WATER QUALITY - ECOLOGICAL MODEL 


TO THE OCONEE RIVER SYSTEM 


By 

Donald J. Smith 
Tetra Tech, Inc. 

3700 Mt. Diablo Boulevard 
Lafayette, California 94549 


BACKGROUND AND PURPOSE 

The U.S. Army Corps of Engineers Hydrologic Engineering Center 
(HEC) is adapting a dynamic water quality-ecological model to a portion 
of the Oconee River in Georgia. While the model is designed to calcu¬ 
late the population dynamics of algae, zooplankton, benthic animals and 
fish, detailed calibration of that portion of the model is beyond the 
limited scope of this project. Since the primary purpose of the model 
will be to evaluate the transient water quality impact of storm runoff 
and waste discharges, the organic sediment and the biological section 
of the model will remain constant during the simulation. 

In lieu of modeling those parameters which were held constant, 
pertinent reports of water quality and biological surveys were reviewed 
to estimate their values. 

This brief report documents the findings of this review. 


RIVER SYSTEM 

The study is limited to the upper reaches of the Oconee River sys¬ 
tem near Athens. Included is the Oconee River between the Barnett 
Shoals Dam and the confluence of the Middle Oconee and North Oconee 
Rivers and approximately twenty-five (25) miles of the Middle Oconee 
River and twenty (20) miles of the North Oconee River. 
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The Middle and North Oconee Rivers are typically 1/2 to 5 feet 
deep and from 50 to 100 feet wide. The Oconee is also typically 50 to 
100 feet wide with depths up to 8 feet. The average gradient is approxi¬ 
mately 4 feet per mile and velocities are characterized as slow to 
moderate. 

Urban development is limited to the Athens area. The remaining 
watershed is rural with many wooded areas. A total of approximately 
seven (7) million gallons a day of municipal and industrial waste water 
is discharged to the Middle Oconee and North Oconee after secondary 
treatment. During periods of high runoff, significant amounts of or¬ 
ganic detritus and sediment are washed into these rivers from the 
watershed. 


WATER QUALITY 

For purposes of characterizing the water quality, the river system 
can be divided into two sections. Section One includes the Oconee River 
and those portions of the Middle Oconee and North Oconee below the two 
Athens sewage treatment plant outfalls. Water quality in this section 
is influenced by the Athens sewage treatment plant effluent. Dissolved 
oxygen is lower and plant nutrients, BOD, and total organic carbon are 
higher than in Section Two. Water quality in Section Two, the remaining 
portion of the study area, is reasonably good. The water of both sec¬ 
tions is quite turbid during periods of high flow. Levels of selected 
water quality parameters reported by state and federal agencies are 
summarized in Table 1. 


ALGAE AND ZOOPLANKTON 

No chlorophyll a data or other direct measurement of suspended al¬ 
gae are available. Some attached algae (periphyton) and the macrophyte 
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Podostemum have been observed where suitable rock substrate is avail- 

p 

able. Nelson (1962) reports Podostemum levels of 10 to 15 g/m (dry 

weight) on the Middle Oconee where ideal substrate conditions exist and 

none on sand or mud substrates. Suitable substrate (bedrock and cobbles) 

should exist where velocities are sufficiently high to prevent deposition 

of sand and silt. If we assume velocities are sufficiently high with 

bottom slopes of 1 foot in 200 feet and average bottom slope of 4 feet 

in one mile, approximately 15% of the substrate is suitable. Fifteen 

percent of the densities reported by Nelson yield average macrophyte 

2 

and periphyton densities from 1.5 to 2.25 g/m . 

An examination of dissolved oxygen data indicates that algal 
photosynthetic oxygen production is not significant. Dissolved oxygen 
never exceeded saturation and no diurnal variation was evident. No de¬ 
pletion of plant nutrients was observed in the data. Nelson also re¬ 
ported that no detectable differences in dissolved oxygen was observed 
between upstream and below his study area. 

For modeling purposes, both algae and zooplankton concentrations 
can be presumed low, near zero. 


BENTHIC ANIMALS 

The make up of the benthic animal (microinvertebrates) population 
has been studied by submerging limestone substrate (LSS) in the water 
for two months. The results of these studies indicate that benthic 
animals can survive throughout the study area if suitable substrates 
exist. Unfortunately the test results do not include total biomass or 
LSS surface areas, therefore, population densities cannot be determined. 

Nelson (1962) observed benthic animal population densities of 2 to 

10 g/m (dry weight) with ideal substrate conditions. If we also assume 

fifteen percent (15%) of the natural substrate is suitable for benthic 

2 

animals, average densities of .3 to 1.5 g/in can be expected. 
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FISH 


Fish population data within the study area are limited to one sam¬ 
pling event in June 1959. Total fish mass was reported. However, the 
length of stream sampled was omitted making it impossible to calculate 
fish biomass per mile. All fish collected were warm water species. 
Approximately fifty percent (50%) were bottom feeding fish. 

Streams of this type typically have a fish population of 100 to 
300 lb/acre (wet weight). Assuming an average channel width of 
75 feet, approximately 600 feet of channel has a surface area of one 
acre. Converting to dry weight per mile, a total biomass of 90 to 260 
lbs/mile or 40 to 120 kg/mile is obtained. 

DETRITUS AND ORGANIC SEDIMENT 

Total organic carbon in the water ranged from 2 to 7 mg/1 C. The 
detritus concentration in the water is generally twice the organic 
carbon level or 4 to 14 mg/1. The detritus level is generally a function 
of flow rate, increased detritus occuring with increased flow. Nelson 
(1962) attempted to correlate river discharge with detritus volatile 
solids) with some success. He typically measured volatile solids of 
1 to 10 mg/1 at moderate flows and up to 50 mg/1 during high flow periods. 

Nelson also measured settleable plant and animal detritus and 
reported typical values of 12 to 20 g/m^. Below the Athens STP outfalls, 
suspended detritus and organic sediment can be expected to increase. A 
10 to 20 percent increase in the above values seems appropriate. 


Nelson, Daniel J., and Scott, Donald C., 1962, Role of detritus in the 
productivity of a rock-outcrop community in a piedmont stream. Depart¬ 
ment of Zoology, University of Georgia. 
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DUALITY PROFILES 
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